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Abstract 

Hall magnetohydrodynamics (MHD) properties near a two-dimensional (2D) X-type mag- 
netic neutral line in the steady state are considered via heuristic and rigorous develop- 
ments. Upon considering the steady-state as the asymptotic limit of the corresponding 
time- dependent problem and using a rigorous development, Hall effects are shown to be able 
to sustain the hyperbolicity of the magnetic field (and hence a more open X-point configura- 
tion) near the neutral line in the steady state. The heuristic development misses this subtle 
connection of the steady state with the corresponding time- dependent problem and predicts 
only an elongated current-sheet configuration (as in resistive MHD). However, the heuristic 
development turns out to be useful in providing insight into the lack of dependence of the 
reconnection rate on the mechanism breaking the frozen-in condition of the magnetic field 
lines. The latter result can be understood in terms of the ability of the ions and electrons 
to transport equal amounts of magnetic flux per unit time out of the reconnection region. 
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1 Introduction 



In resistive magnetohydrodynamics (MHD) the ion inflow is the only means to transport 
magnetic flux into the reconnection layer (Knoll and Chacon[l]). As the resistivity is de- 
creased, large magnetic pressure gradients which develop upstream of the reconnection layer 
start inhibiting the ion inflow and the magnetic flux transport into the reconnection layer 
(and hence the reconnection rate) - the so-called pressure problem (Clark [2]). The Hall 
effect (Sonnerup [3]) can overcome the pressure problem (Dorelli and Birn [4], Knoll and 
Chacon [5] ) , thanks to the decoupling of electrons from ions on length scales below the ion 
skin depth rfj. So, if the reconnection layer width is less than rfj, the electron inflow can keep 
on going which transports the magnetic flux into the reconnection layer and hence reduces 
the flux pile-up. Previous numerical work (Shay et al. [6], Rogers et al. [7], Knoll and 
Chacon [5]) indicated that the dissipation region in Hall MHD, as increases, changes from 
an elongated current sheet geometry (Sweet [8]-Parker [9] type) to a more open X-point 
geometry (Petschek [10] type). However, recent fully kinetic simulations (Daughton et al. 
[11], Karimabadi et al. [12]) and EMHD-based treatments (Chacon et al. [13]) have shown 
that elongated current sheets are also possible. On the other hand, more recent particle-in- 
cell (PIC) simulations (Shay et al.[14]) show spatial localization of the out-of-plane current 
to within a few d'-s of the X-line, so this controversy continues. It is therefore in order 
to shed further light on this issue. In this paper, we consider Hall MHD properties near 
a two-dimensional (2D) X-type magnetic neutral line in the steady state via heuristic and 
rigorous developments and investigate whether or not the Hall effects favor the hyperbohcity 
of the magnetic fleld near the neutral line. The heuristic development turns out to be useful 
in providing insight into the lack of dependence of the reconnection rate on the mechanism 
breaking the frozen-in condition of the magnetic fleld lines. 



2 Governing Equations for Hall MHD 



Consider an incompressible, two-fluid, quasi-neutral plasma. The governing equations for 
this plasma dynamics are (in usual notation) - 



nrrii 



dt 

'dvj 
dt 



+ (Ve • V)V, 

+ (vi • V)vi 



— Vpe ~ neCE -\ — Vg x B) -|- ner^J 
c 

— Vpi -|- ne(E H — Vj x B) — nerjJ 



(1) 

(2) 



where. 



V • V, = 

V • Vi = 
V-B = 

V X B = -J 

c 

^ ^ l^B 

VxE = — — 
c at 

J = ne{vi - Ve). 



(3) 
(4) 
(5) 

(6) 
(7) 

(8) 
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Neglecting electron inertia (me — > 0), equations (1) and (2) can be combined to give an 
ion equation of motion - 



nrrii 



^ + (v.-V)v. 



-V(pi+Pe) + -JxB (9) 
c 



and a generalized Ohm's law 



E + -Vi X B = 77J + — J X B. (10) 
c nec 

Non-dimensionalize distance with respect to a typical length scale a, magnetic field with 
respect to a typical magnetic field strength Bq, time with respect to the reference Alfvcn 
time ta = o/Vaq where Vaq = Bo/y/rriin, and introduce the magnetic and velocity stream 
functions according to 



B = Vip X + 6i 
Vi = V0 X + w'l 



z 



(11) 



and assume the physical quantities of interest have no variation along the 2;-direction. The 
Hall magnetic field b is believed to be produced by the dragging of the in-plane magnetic 
field in the out-of- plane direction by the electron near the X-type magnetic neutral line ([6], 
[7]). Equations (9) and (10), then yield 

^ + [^>] + a[6,V']=r?VV (12) 

db 

- + [b, 0] + a[ij, VV] + w] = fjW'b (13) 

|(V20) + [V20,0] = [V^V',V'] (14) 

^ + = (15) 



where, 



[A, B] = VAxVB ■ i, 

= — , V = — —■ 
a 



3 A Heuristic Analysis 

It is instructive to do a heuristic analysis to develop an estimate on the geometry of the 
dissipation region prior to a more rigorous formulation. Let the dissipation region have a 
length L in the outfiow x-direction and a width 6 in the infiow |/-direction. 
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We then have from equation (5), 



Bx By 



L S 

We have from the ;2-component of equation (10), 

VJz ~ JxBy 

nec 

and on using equation (6), (17) becomes 



or 



Bx I B-^ 

Vc—. ^^^y 

nec 



rjBx —ByB^. 

nec 



Next, the z-component of the curl of equation (10) gives 

77cV'S,= — (B-V)J, 

nec 

from which, 

Bz 1 Bx Bx 



TIC- 



S'^ nec L S 



and on using (16), (20) becomes 
or 

(18) and (21) give 
Using (22), (18) gives 

Using (16), (23) leads to 
On rearranging, (24) gives 
where. 



1 By Bx 



S'^ nec S S 



vBz —BxBy. 

nec 



Bz ~ Bx. 



By ~ necTj. 

j^Bx^ necrj. 

6 1 
L aS 



(25) was also given by Chacon et al. [13] by using a more rigorous formulation. 
Noting that in the Hall resistive regime 



(7 > 1, S»l 

(25) imphes that the diffusion region in the Hall resistive region may be expected to be elon- 
gated. Rigorous formulation in Section 4, on the other hand, shows that this result represents 
only part of the story because the heuristic development fails to recognize the steady state as 
the asymptotic limit of the corresponding time- dependent problem (Shivamoggi [15]). The 
latter aspect makes the actual story more complicated than what is conveyed by the heuristic 
development. 

It is interesting to note, however, that the above heuristic analysis sheds some light on the 
conjecture (Mandt et.al. [16], Shay and Drake [17]) that the reconnection rate in Hall MHD 
is primarily controlled by ions (which are decoupled from the electrons) and is independent 
of the mechanism that breaks the frozen-in condition of the magnetic field lines (resistivity 
or electron inertia). 

For the Hall-resistive case, the reconnection rate is give by 



E ~ r/ J, (26) 
and on using equation (6) and (23), (26) becomes 

E ^ Ac:^. (27) 
nec 6 

Using (16), (27) becomes 

E ^. (28) 

neL 

On the other hand, if we consider the electron inertia to constitute the mechanism that 
breaks the frozen-in condition of the magnetic field lines, the Ohm's law now takes the form 
(Coppi et al. [18]) 

E+-VeXB = — — . 29 
c uj± at 

Pe 

In the electron-inertia case, the reconnection rate is therefore given by 

E (30) 

a;„2 at 

Using equation (6), (30) may be rewritten as 

1 1 5x 

E ^. (31) 

UJp2 Ta^ 



where, 

L _ B 



TAe ~ T^, = 



X 



Va, ' ^/nme 
Taking 5 ~ de, (31) becomes 



5 



Er^—^ (32) 

neL 

which is the same as the one, namely, (28), for the Hall resistive case! This appears to 
support the conjecture (Mandt et al. [16], Shay and Drake [17]) that the reconnection rate 

is independent of the mechanism that breaks the frozen-in condition of the magnetic field 
lines. A similar conclusion was reached by Chacon et al [13] who considered the electron 
hyperresistivity to constitute another mechanism that breaks the frozen-in condition of the 
magnetic field lines. 

It is of interest to note that (28) and (32) may be rewritten as 

E^J^^bJY-\(^\^bJ'M(^-1]. (33) 



neL ^V*^ J \^ J ^V*^ J \^ , 

(33) shows that the lack of dependence of the reconnection rate on the mechanism breaking 
the frozen-in condition of the magnetic field lines can be understood in terms of the ability 
of the ions and electrons to transport equal amounts of magnetic flux per unit time out of 
the reconnection region. The insensitivity of the reconnection rate, according to (33), on 
the particle mass has been confirmed by the recent partial-in-cell simulations ([14]). We are 
assuming here, as confirmed by the recent numerical simulations (Drake et al. [19]) that, in 
the electron-inertia case, the outflow velocity of the electrons from the dissipation region is 
given by the Alfven speed based on the upstream magnetic held with the width of the 
dissipation region scaling with d^. 

4 Steady-state Properties Near an X-type Neutral Line 

Consider Hall MHD properties near a 2D X-type magnetic neutral line in the steady state. 
Equations (12)-(15) now become 

-cE + ^10001 - ^01010 + (y{hioipm - &oi^io) = fj{ip2o + 1P02) (34) 

^'lO^Ol - + O-['01o('021 + -003) - V'0l(V'30 + "012)] + i^wW^i - ^'OI^^IO = fjiho + ^'02) (35) 

01o(021 + 0O3)-0Ol(03O + 012)-[^/'lo(^21+^O3)+^Ol(^3O + ^12)] = -^^(^40 + 2022 + ^04) (36) 
W1O0O1 - WO101O + V'lO&Ol - V'Ol&lO = 1^{W20 + W02) (37) 



where. 



Qm+np I 



dx™-dy^ c dt 

and we have now included in equations (14) and (15) viscous effects in the plasma which 
become important near the magnetic neutral line (Tsuda and Ugai [20]); u is the viscosity 
coefficient. 
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Following Cowley [21] and Shivamoggi [22], let us expand the velocity and magnetic fields 
in a Taylor series about the neutral line taken to be at a; = 0, y = 0. Equations (34)-(37) 
may then be used to derive relationships between the coefficents of the series. The latter 
are simply the partial derivatives of the velocity and magnetic fields at the neutral line. 
Motivated by the symmetry properties of equations (34) - (37) in the ideal limit, we may 
consider i/j and w to be even functions of both x and y, and and b to be odd functions of 
both X and y (this also enables the out-of-plane magnetic field b to exhibit the quadrupolar 
structure (Terasawa [23]) characteristic of Hall MHD - this has also been confirmed by 
laboratory experiments (Ren et al. [24]). Thus, we write 



2n" 



(2m)!(2n)! 



= E E <^2^+W^2m+l)!(2n + l)! 

m n 

„2m+l„,2n+l 



x^-'-y 

52m+l,2n+l 



(2m + l)!(2n + l)! 

j.2my2n 

(38) reflects the fact that the origin in the x, y — plane is both the X-type neutral point and 
a stagnation point of the flow. 

Using (38), equations (34) and (37) give, on evaluation at the origin, 

r?(^2o + ^02) = -cE (39) 

iy{W2o + W02) = (40) 

Let us differentiate equations (34) and (37) with respect to x and y separately, and 
differentiate the resulting four equations, respectively, with respect to x, y and x, y. We then 
obtain 

'030001 + 2-020011 + "010021 " "021010 " 2-011020 " "001030 

+ a{b3Qipoi + 262oV'ii + ^'io'02i - &2iV'io - 26iiV'2o - feoiV'3o) 

= 77(V'40 + -022) (41) 
■012001 + 2-011002 + -010003 " "003010 " 2-002011 " -001012 
+ Cr{bi2lpoi + 26llV'02 + &loV'03 - &03V'lO - 26o2V'll - ^'OlV'12) 

= r/(^22 + ^04) (42) 

W3O0O1 + 2^20011 + ^^10021 - W2l(l)l0 - 2^11020 - ^01030 
+ 0'3O&O1 + 20^02^11 + ^10^21 - ^21^10 - 2^^11620 - ^^01^30 

= Z^(W40 + W22) (43) 
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+ 2^11002 + W1O0O3 - ^t!O301O " 2^02011 - ^01012 
+ Ipubol + 2-011602 + -010^*03 - '003^*10 - 2-002^11 - ipoibu 

^ 1'{W22+W04). (44) 

Next, let us differentiate both equations (35) and (36) with respect to x and then both 
with respect to y. We then obtain 

&210O1 + &2O0O2 + &1O012 " &1201O " &O202O " &O1021 + 0-[^2l(^21 + ^03) + ^2o(^22 + ^04) 

+ 0^10(0^32 + 0^14) - 0^12(030 + -012) - 0^02(0^40 + 022) - 0^01(041 + 0'23)] 
+ 0^21^01 + -020^02 + 0^10^12 " -012^^10 " 0^02^20 - 0^01^21 = r)(&31 + &I3) (45) 
021 (021 + 003) + 020 (022 + 004) + 0io(032 + 014) " 012 (030 + 012) " 002 (040 + 022) " 001 (041 + 023) 
= 0'21 (0'21+0'O3)+V'2O (0'22+0'O4)+0'lO (0'32+0'l4) -0'l2 (0'3O+V'l2) -0'O2 (0'4O+V'22) "V'Ol (0'41+0'23) 

-1/(051 +2033 + 015). (46) 

We now use (38) and evaluate equations (41) - (46) at the origin: 

?7(*4o + *22) = 2*2o(^ii - crBn) (47) 

77(*22 + *04) = -2*02(^11 - crBu) (48) 

77(^31 + S13) = (7[*20(*22 + *04) - *02(*40 + ^22)] 

+ *20W"02 - ^02^^20 (49) 

*20(*22 + *04) - ^02(^40 + ^22) = i^($51 + 2$33 + $15) (50) 

1/(^40 + W22) = 2W^20*11 + 2^02511- (51) 

iy{W22 + W04) = -21^02^11 - 2^20^11 (52) 
Using equations (47) and (48), equation (50) gives 

4($ii - (7Sii)*20*02 = -i^^($51 + 2$33 + ^is)- (53) 

Equation (53) shows that, in the MHD resistive viscous case ((7 = 0, fj and u 7^ 0), one 
has 

^20 ^ 0, ^02 7^ (54) 



8 



so the magnetic field to lowest order can be hyperbolic (Shivamoggi [22]). On the other hand, 
in the inviscid or non-resistive case, when Hall effects are included {Bu 7^ - Hall effects 
materialize only via their signature - the quadrupolar out-of-plane magnetic field pattern), 
on considering the steady state as the asymptotic limit of the corresponding time- dependent 
problem^ (Shivamoggi [15]), we have from the latter formulation 

$11 - aBii = 0. (55) 

Using (55), equation (53) shows that, even in the inviscid limit {u ^ 0), equation (54), thanks 
to Hall effects, continues to be valid. So, Hall effects are able to sustain the hyperbolicity 
of the magnetic field (and hence a more open X-point configuration) near the neutral line. 
Equation (55) implies that the level curves of the out-of-plane magnetic field are also the 
streamlines of the in-plane ion flow. 

Further, using (55), equations (49), (51) and (52) lead to the following compatibility 
conditions on the Taylor expansion coefficients of the out-of-plane components of the velocity 
and magnetic fields - 

r)(fi3i + 5i3) = (56) 
^02 + crWo2 = (57) 
^20 + (yW2o = 0. (58) 

5 Discussion 

In this paper, we have considered Hall MHD properties near a 2D X-type magnetic neutral 
line in the steady state via heuristic as well as rigorous developments. Upon considering the 
steady state as the asymptotic hmit of the corresponding time- dependent problem and using 
a rigorous development. Hall effects are shown to be able to sustain the hyperbolicity of the 
magnetic field (and hence a more open X-point configuration) near the neutral line in the 
steady state. The heuristic development misses this subtle connection of the steady state 
with the corresponding time- dependent problem and predicts only an elongated current-sheet 
configuration (as in resistive MHD). On the other hand, this development also shows that 
the electron-hyperresistivity effects are similar to that of ion viscosity p and allow for the 
possibility of a X-point magnetic field configuration (see Appendix) as found also by Chacon 
et al. [13]. The heuristic development, however, turns out to be useful in providing insight 
into the lack of dependence of the reconnection rate on the mechanism breaking the frozen-in 
condition of the magnetic field lines. The latter result can be understood in terms of the 



^Indeed, one way to resolve difficulties that arise in insuring uniqueness of solution of 
stationary linear wave problems is to pose a more realistic unsteady problem with suitable 
initial conditions applied at some finite time t = —to, say, in the past, and then letting 
to ^ 00 in the solution (Lighthill [25]). 
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ability of the ions and electrons to transport equal amounts of magnetic flux per unit time 
out of the reconnection region. 
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7 Appendix 

On including the effects of ion viscosity and electron hyperresistivity, equations (12)- (15) 
become 

^ + [V', 4>] + a[h, i^] = ?7V V + ?7ft VV (59) 

dh 

- + [6, 4>] + a[^, VV] + [V', w] = fiV% + r]hV% (60) 
^( W) + [VV, 0] + [V', VV] = i/VV (61) 

^ + [w, 0] + [V', h\ = vV'w (62) 

where r]h is the electron hyperresistivity coefficient. 

Following through the development outlined in Section 4, we have in place of equations 
(41), (42) and (45), respectively, 

^30001 + 2^20011 + ^10021 - ^/'2101O " 2^/^ii02O " ^01030 
+ (7-(630^01 + 2620^11 + &10^21 - &21^10 - 26ii^/'20 " &01^3o) 

= ?7(V'40 + V'22)+?7/i(V'60 + V'24) (63) 

V'12001 + 2V'110O2 + V'10003 - V'03010 " 2V'O2011 " V'O1012 
+ O'(fcl2V'01 + 26iiV'02 + &loV'03 - ^'OsV'lO - 26o2V'll - ^'OlV'12) 

= fi{i^22 + V'04) + ?7ft(V'42 + V'oe) (64) 

&210O1 + &2O0O2 + &10</'12 - ^12010 " &O202O - &O1021 + Cr[^2l(^21 + ^03) + ^/'2o(V'22 + ^"04) 



10 



+ ^^10(^^32 + Ipu) - ^^12(^^30 + Ipu) - 1p02{lp40 + 1P22) - V'0l(V'41 + ^^23)] 



+ V'2i«^oi + V'2o«^02 + ipwWi2 - ipuwio " ^'02^^20 " V'oi«^2i = vihi + ha) + Vhihi + hb)- (65) 
Using (38) and evaluating equations (64) - (66) at the origin, 

vi^Ao + ^22) + Vhi^GO + ^24) = 2^20 ($11 - (^Bu) (66) 

r}(*22 + *04) + Vhi'^42 + ^oe) = -2*02(^11 - c^^ii) (67) 

77(^31 + Bi3)+r]h{B5i + S15) = (7[*20(*22 + *04) - *02(*40 + ^22)] + ^20^^02 - *02W"20- (68) 

Using equations (67) and (68), equation (50) gives 

4($n - (7Sn)*2o*02 = -i^f}{^5i + 2^33 + $15) - ?7h(*60 + *42 + *24 + ^og)- (69) 

It is of interest to note that the right hand side in equation (70) comes from the term - 



+ 



^ dx^ ^ dxdy^ ) \dx'^ dy'^ 



which highhghts the similarity between the ion-viscosity and electron hyperresistivity con- 
tributions - both allow the possibility of a X-point magnetic field configuration. 
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